The diffusion coefficient of propagating fronts with multiplicative noise 
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Recent studies have shown that in the presence of noise both fronts propagating into a metastable 
state and so-called pushed fronts propagating into an unstable state, exhibit diffusive wandering 
about the average position. In this paper we derive an expression for the effective diffusion coefficient 
of such fronts, which was motivated before on the basis of a multiple scale ansatz. Our systematic 
derivation is based on the decomposition of the fluctuating front into a suitably positioned average 
profile plus fluctuating eigenmodes of the stability operator. While the fluctuations of the front 
position in this particular decomposition are a Wiener process on all time scales, the fluctuations 
about the time averaged front profile relax exponentially. 

PACS numbers: 05.40.-a, 47.54.+r, 05.45.-a 



I. INTRODUCTION 

One of the aspects of front propagation that have 
been studied in the literature in recent years is the 
effect of fluctuations on propagating fronts Q . In 
particular it has been found that in the presence of 
noise both one-dimensional fronts between a stable and a 
metastable state ( "bistable fronts" ) and so-called pushed 
fronts which propagate into an unstable state |^] , exhibit 
a diffusive wandering about their average position [Q. 
This contrasts with the fluctuation behavior of so-called 
pulled fronts propagating into an unstable state which is 
subdiffusive ||. In this paper we shall consider only the 
case of pushed and bistable fronts, however. 

Recently, Armero et al. j| derived an expression for 
the effective diffusion coefficient of a pushed front in the 
stochastic field equation 



dt dx 2 

with a noise term whose correlations are 



f](x,t) = 0, 



ri(x,t)r)(x',t') = 2C(\x - x'\/A)S(t - t'). 



(1) 



(2) 



(3) 



In (jlj) , / is a nonlinear function of the field (f) with a sta- 
ble state at <fi — 1 and either a (meta)stable or unstable 
state at <p — and g{4>) is some other general nonlinear 
function. In (||) and (|^) the overbar denotes an average 
over the realizations of the noise. In order that our noise 
of Stratonovich type is well-defined, we have introduced 
a spatial cutoff in the noise correlation function (^) (see 
[pj for further details). 

The derivation in Q] of the effective front diffusion coef- 
ficient D f relied on a small-noise stochastic multiple scale 



analysis that was based on the idea that the mean square 
displacement of the front about its average position was 
slow relative to the deterministic relaxation of the front. 
The basic idea was that only the low-frequency compo- 
nents of the noise are responsible for the front wandering, 
so that the high-frequency components, which renormal- 
ize the front shape and its velocity, could be implicitly in- 
tegrated out. This led to an ansatz for the relative scaling 
of fast and slow time variables where the small parame- 
ter governing the separation of time scales was the diffu- 
sion coefficient Df of the front itself. The method then 
sclfconsistently provided an explicit prediction for Df. 
The main weakness of the approach was that the above 
coarse-graining procedure could not be carried out ex- 
plicitly, since while there is a separation of time scales for 
the average quantities, a scale separation scheme is not 
natural for the fluctuating quantities. Hence the deriva- 
tion had to rely on an uncontrolled ansatz. In this Brief 
Report we therefore reconsider this problem. We justify 
the previously derived result for Df with a systematic 
small-noise expansion based on decomposing the motion 
of the front into a diffusive motion of the properly de- 
fined front position plus fluctuations about the average 
front profile. Technically the fluctuating front position 
is defined by requiring that the fluctuations about the 
mean front profile are orthogonal to the (left) translation 
mode. This derivation shows that the previous multiple 
scale ansatz is not quite adequate, and it will clarify the 
connection between the separation of time scales invoked 
in Ref. [Q, the small noise expansion and the existence 
of a finite gap in the linearized evolution operator. The 
key point of our new derivation is the fact that there is 
a unique choice for the collective coordinate X{t) of the 
front profile to be a memoryless Markovian process, and 
that the fluctuations about the average profile then relax 
exponentially. This relaxation can be deduced from the 



f 



spectrum of the linearization operator about the average 
front profile. In addition our method provides a general 
strategy to address the problem of fluctuations of fronts 
and other coherent structures. 



II. DERIVATION OF THE EFFECTIVE 
DIFFUSION COEFFICIENT 

We can rewrite Eq. ([!]) in terms of a noise term R whose 
average R is zero and a deterministic renormalized part, 



i) 2 



h((j,) + s 1/2 R(4>,x,t), 



(4) 



dt dx 2 

using Novikov's Theorem, as discussed in jij. In Eq. 

h((f>)=f(cl>)+eC(0)g'(<j>)g(<t>), (5) 
R{4>, x, t) = g(<f>)r}(x, t) - e 1/2 C(0) 5 ' '(<P)g(<t>), (6) 

where C(0) is of order A -1 , so that Eq. (||) yields a 
delta correlation in space in the limit A — » 0. The 
main idea of the derivation is to introduce a collective 
coordinate X(t) for the position of the front. Of course 
there are various choices for the position X(t), but as we 
shall show a particular choice makes the equations quite 
transparent. We decompose the fluctuating field cj> as 



= Mt-x(t)) + MZ-x(t),t). 



(7) 



Here <fio is the solution of the ODE for the shape of a 
deterministic front with velocity vr, the velocity of the 
deterministic front associated with Eq. (^) with R = 
(the subscript R on vr reminds us that the front speed is 
determined by h(<fi) rather than /(</>), and hence is renor- 
malized due to the presence of the noise). In other words, 
4>q satisfies 



VR- 



di 



(8) 



While 4> is a non-fluctuating quantity, <j>i is a stochastic 
field which contains the fluctuations about (j)Q. In the 
above, £ = x — VRt is the proper variable for a deter- 
ministic front moving with the asymptotic velocity vr, 
but note that in (Q) the fields are written in terms of the 
shifted variable 



6f =£-X(t) =x-v R t + X(t), 



(9) 



where X(t) is the rapidly fluctuating front position whose 
explicit definition in terms of a spatially averaged front 
profile is given below. 

As is well known, the derivation of a moving boundary 
approximation for deterministic equations (see e.g. [^|J^] 
and references therein) normally proceeds by projecting 
onto the zero mode. Indeed associated with the front 
solution 0o of (p|) is a zero mode of the stability operator 



C 



vr-q^ + h'{cj>o), 



(10) 



which is obtained by linearizing about 4>q. This zero 
mode expresses translational invariance, and indeed im- 
plies that 



(o) 







(0) 



#0 



(11) 



In our case the operator C is not self-adjoint, since 
vr ^ 0; as a result the left eigenmode $^ is different 
from ^r , but it is known to be (see e.g. |4|||): 



£+4 0) = 



(o) _ 



If- (12) 



As we mentioned above, a particular definition of the 
position X{t) is especially convenient: we take X(t) de- 
fined implicitly by the requirement that the fluctuating 
field </>i is orthogonal to the left zero mode. Indeed, defin- 



ing 



/oc 
d£A(0B(0, 
-oo 



(13) 



we require 



<<fi 0) <M^)> 

= J d^^{<t> - - X{t)) = 0- (14) 

Note that at any moment, the fluctuating front position 
X(t) is defined in terms of weighted spatial average of the 
fluctuating field (f>. 

Upon substitution of Eq. (^) into (^) and linearization 
in <pi (which is justified for small noise), we obtain 

^ = £4x-X{t)^-+R{ f h,i i t). (15) 
dt de,x 

Note that we have also approximated R((f>, £,x,t) by 
R((j>Q^,t), which again is correct to lowest order in the 
noise. 

In addition to the zero mode, the operator C will in 
general have right eigenmodes with eigenvalues —oi 

= -a,*® i^q (16) 

and with associated left eigenfunctions $^ = e VR ^^>^. 
Our convention to have the eigenvalues — ai anticipates 
that the dynamically relevant front solution is stable, 
so that all eigenvalues 07 are positive. Moreover both 
for fronts between a stable and a metastable state and 
for pushed fronts propagating into an unstable state, the 
spectrum is known to be gapped [ [l0|]Tl| ] , i. e. the smallest 
eigenvalue is strictly greater than zero [[10|Jll| . 
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Since <f>\ is orthogonal to we can expand cf>i in 

(0 



terms of the eigenmodes 3> R (I > 1) of C as 



(17) 



Substitution of this expansion into (|l^) then yields upon 



projection onto the zero mode $ 



(0) 



(18) 



Taking the square of this result, integrating and averag- 
ing over the noise, 



X 2 (t) = 2D ft = dt' dt" X(t')X(t"), (19) 
J a Jo 



then yields with (g), @ and @ 

J^ e 2 ^0oM)V(0o) 



D, = e* 



[/de e «««(d^o/de) 



212 



(20) 



This is precisely the result derived earlier in [|4j. To 
lowest order in the present small-noise expansion, the av- 
erage front profile is simply <pQ. However, notice that <po 
contains a dependence on A through 0(0) in h{<j>). The 
parameter C(0) must be considered as an independent 
one, so that the result (^0|) has to be interpreted as to 
first order in e but to all orders in e/A. 

The above derivation allows us to also obtain the re- 
laxation of a fluctuation about the average. Indeed, 
upon substituting (|l^) into fll5| ) and projecting onto the 
left zero modes, using {^^Q^} = S nrn for normalized 
eigenmodes, we obtain to lowest order 



dai 
dt 



-criai 



(21) 



as terms X (t)d<fii / d£ are of higher order in e. Note that 
each mode is damped and has its noise strength weighted 
by • One can derive from here in a straightforward 
way the mean square fluctuations about the average pro- 
file. 

We finally note that our discussion clarifies the diffi- 
culty of using a separation of time scales argument for 
the derivation of the effective diffusion coefficient: the 
collective coordinate X(t) is a memoryless Markov pro- 
cess, and hence the changes in the position have zero cor- 
relation time while the average of X 2 (t) changes slowly. 
The coefficients ai(t), on the other hand, have a finite 
correlation time and hence are correlated on timescales 
in between the one of instantaneous position X{t) and 
the mean square wandering X 2 (t). 



III. CONCLUDING REMARKS 



We have reported an improved derivation of the diffu- 
sion coefficient of propagating pushed fronts with multi- 
plicative noise, previously found in Ref. The present 
derivation is more transparent and elegant since it is fully 
explicit and based on standard projection techniques. 
The key point is the identification of a definition of the 
front position which naturally implies the diffusive wan- 
dering of the front, and avoids invoking an uncontrolled 
hypothesis in addition to the basic assumption of small 
noise strength. This has also clarified that the time scale 
separation used in Ref. can be traced back to the 
small noise approximation together with the existence of 
a finite gap in the spectrum of the linearized evolution 
operator. All these considerations can be generalized to 
the effect of fluctuations on other types of coherent struc- 
tures. 

Our derivation of the solvability expression (|2^) for 
Df of a propagating front shows that the collective co- 
ordinate X(t) responds instantaneously to the noise R: 
There is no memory term in (|l|), so that X(t) is Marko- 
vian and, more precisely, it coincides with the Wiener 
process (to lowest order in the noise strength) . We stress 
that this is only true for our particular definition of X(t) 
in terms of the orthogonality of 4>i to the left zero mode. 
For any other definition, like the usual one to define the 
front position as X(t) — j d£ </>(£), X(t) will not be a 
Markov process, and would show only diffusive behavior 
at sufficiently long time scales. 

As a byproduct of our derivation we have also ob- 
tained an explicit expression for the relaxation behavior 
of the fluctuations about the mean front profile. Not sur- 
prisingly, the larger the gap in the spectrum, the faster 
the relaxation. As is well known, in models in which 
there is a transition from the pushed regime to the pulled 
regime, the gap closes upon approaching the transition 
from the pushed side . Hence the relaxation becomes 
slower and slower. As is discussed in 10 , in the pulled 
regime the spectrum is gapless and this leads to anoma- 
lous power law relaxation of deterministic fronts towards 
their asymptotic speed and shape. As a result pulled 
fronts cannot be described by a moving boundary ap- 
proximation H and in the presence of fluctuations they 
exhibit subdiffusive wandering |jj in one dimension and 
anomalous scaling in higher dimensions [|l2|,[l3| . 
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